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Abstract 

O-J We consider a tree-level enhancement to the violation of lepton flavour universality in 

light meson decays arising from modified W tv couplings in the standard model minimally 
QJ , extended by sterile neutrinos. Due to the presence of additional mixings between the active 

(left-handed) neutrinos and the new sterile states, the deviation from unitarity of the leptonic 
mixing matrix intervening in charged currents might lead to a tree- level enhancement of Rp = 
£SJ ■ r(P — > ev)/T(P —¥ (J,v), with P = K, ir. Wc illustrate these enhancements in the case of the 

inverse seesaw model, showing that one can saturate the current experimental bounds on Arx 
(and Arv), while in agreement with the different experimental and observational constraints. 
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1 Introduction 



Lepton flavour universality (LFU) is one of the distinctive features of the Standard Model of 
strong and electroweak interactions (SM); hence, any deviation from the expected SM theoretical 
estimates in electroweak precision tests will signal the presence of New Physics (NP). Here we 
focus on light meson (K and it) leptonic decays which, in view of the expected experimental 
precision, have a unique potential to probe deviations from the SM regarding lepton universality. 

In the SM, the dominant contribution to T(P — > tv) (P = K, ir) arises from W boson mediated 
exchanges. The prediction of each specific decay is heavily plagued by hadronic matrix element 
uncertainties; however, by considering the ratios 

= T(K+ -> e+v) r(vr+ ->■ e+v) 

K " T(K+ -> fi+v) ' " ~ r(vr+ -> (i+v) ' { ' 

the hadronic uncertainties cancel out to a good approximation, so that the SM predictions can 
be computed with a high precision. In order to compare the experimental bounds with the SM 
predictions, it proves convenient to introduce a quantity, Arp, which parametrizes deviations from 
the SM expectations, possibly arising from NP contributions: 

i? cxp 

R c * p = R s p M (1 + Ar P ) or equivalently Ar P = - 1 . (2) 
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The comparison of theoretical analyses [H[2] with the recent measurements from the NA62 col- 
laboration [3] and with the existing measurements on pion leptonic decays [1] 

i?! M = (2.477 ± 0.001) x 10" 5 , R e x p = (2.488 ± 0.010) x 10" 5 , (3) 

R SM = (1.2354 ±0.0002) x 10~ 4 , Rf p = (1.230 ± 0.004) x 10" 4 (4) 

suggest that observation agrees at lcr level with the SM's predictions for 

Ar K = (4 ± 4) x 10~ 3 , Ar n = (-4 ± 3) x 10~ 3 . (5) 

The current experimental uncertainty in Atk (of around 0.4%) should be further reduced in the 
near future, as one expects to have SRk/Rk ~ 0.1% [5], which can translate into measuring 
deviations Ar^ ~ 0(1O -3 ). There are also plans for a more precise determination of Ar w [6l[7]. 

Whether or not Rp can probe a model of NP naturally depends on the nature and on the 
expected size of the corresponding contributions to Arp. LFU in light meson decays can be 
violated due to (i) a new Lorentz structure in the four-fermion interaction (arising from the 
exchange of new fields); (ii) corrections to the SM charged current interaction Wlv vertex. The 
first possibility has been extensively discussed in the literature, especially in the framework of 
models with an enlarged Higgs sector, since in the presence of charged scalar Higgs, new tree- 
level contributions are expected. However, as in the case of most of Two Higgs Doublet Models 
(2HDM)0, or supersymmetric (SUSY) extensions of the SM, these new tree-level corrections are 
lepton universal [9]. In SUSY models, higher order non-holomorphic couplings can indeed provide 
new contributions to Rp |10H14| . but in view of current experimental bounds (collider, 5-physics 
and r-lepton decays), one can have at most Atk < 10~ 3 in the framework of unconstrained 
minimal SUSY models [13] • Corrections to the Wtv vertex (case ii) can also induce violation of 
LFU in charged currents. Loop level corrections to the latter vertex have been considered and, as 
referred to in |llj . new particles (at a scale ~ Anp) can indeed account for such a contribution, 
but the effect is of order (a/47r) x (m^/Aj^p) and generally well below experimental sensitivity. 

The tree-level corrections to charged current interactions, once neutrino oscillations are in- 
corporated into the SM, are more interesting. In this case, and working in the basis where the 
charged lepton mass matrix is diagonal, the flavour-conserving term oc glj^P^v^W^ now reads 

- c cc = JLui%rPmW^ + c.c. , (6) 

where Ul l is a generic leptonic mixing matrix, i = 1, . . . ,n u denoting the physical neutrino states 
(not necessarily corresponding to the three left-handed SM states = v£) and j = 1,...,3 the 
charged lepton flavour. In the case of three neutrino generations, Ul corresponds to the unitary 
PMNS matrix and flavour universality is preserved in meson decays: since one cannot tag the 
flavour of the final state neutrino (missing energy), the meson decay amplitude is proportional to 
{U v ut)jj = 1, and thus no new contribution to Rp is expected. 

In order to account for neutrino masses and mixings, the SM can be extended with new neutral 
sterile fermion states (n u > 3). In many neutrino mass models there are additional singlet states; 
the existence of sterile states is also strongly supported by current data from reactor experiments, 
cosmology, as well as indications from large structure formation |15| . 

In the presence of sterile states, the Wtv vertex is proportional to a rectangular 3x n u matrix 
Uv, and the mixing between the left-handed leptons vlAl corresponds to a 3 x 3 block of f/j, , 

t^PMNS — > ^PMNS = (1 — V) ^PMNS • (7) 



As recently shown in [S] , the conclusions can be different in a more generalised 2HDM framework. 
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The larger the mixing between the active (left-handed) neutrinos and the new states, the more 
pronounced the deviations from unitarity of C/pmns > parametrized by the matrix 77 [16] . The active- 
sterile mixings and the departure from unitarity of C/pmns can be at the source of the violation 
of LFU in different neutrino mass models which introduce sterile fermionic states (gauge singlets, 
not necessarily right-handed) to generate non-zero masses and mixings for the light neutrinos. If 
these new states have very small mixings with the left-handed neutrinos, then £/pmns ~ C^pmns 
and no effect is expected regarding LFU violation (for instance in the case of high-scale type 
I and III fermionic seesaws). On the other hand, if the model is such that the singlet states 
are not excessively heavy and have large mixings with the active ones, then the deviations from 
unitarity (i.e. rj) can indeed lead to a (tree-level) correction of the Wtv vertex. The impact of the 
non-unitarity of the lepton mixing matrix on leptonic light meson decays was first investigated 
in [17] . prior to the confirmation of neutrino oscillations. In this work, we revisit this idea in the 
light of recent neutrino data and in view of the present (and future) experimental sensitivities to 



Corrections to the Wlv vertex can arise in several scenarios with additional (light) singlet 
states, as is the case of z^SM |18| . low-scale type-I seesaw [19] and the Inverse Seesaw (ISS) [20J, 
among other possibilities. As we proceed to discuss, should the masses of the new singlet states 
(zAi) be light enough (lighter than the decaying meson), or if the large mixings induce significant 
unitarity violation^! (even for m Us 3> mp), then one might have sizeable Arp, within reach of 
current experimental sensitivities. 

In the following section [2] we compute in a model-independent approach Arp in the presence 
of additional fermionic sterile states; we then briefly review in Section [3] the most important ex- 
perimental and observational constraints on the mass of the additional singlet states. In Section 
we consider the case of the inverse seesaw model to illustrate the impact of sterile neutrinos on 
Arp. Our concluding remarks are summarised in Section [5l 

2 Atk in the presence of sterile neutrinos 

Let us consider the SM extended by N s additional sterile states, and conduct a general formulation 
of leptonic light meson decays. The matrix element for the meson decay P — > IjVi has the generic 
form 



no sum implied over the indices of the outgoing leptons i, j. Notice that now one has i = 1, . . . , 3+ 
N s . The expressions for A and B can be read from the effective hamiltonian (obtained after 
integrating out the W boson in Eq. ([6])), 



Arp 



M ij = u l/i (A»P R + B l ^P L )v lj , 



(8) 




(9) 



implying that in this framework one has 



{AT 
(B) ij 



(A 
{B l 



W\ij 




4G F V&MfpUP*m h ; 



(10) 




where fp denotes the meson decay constant and uii- Ui the mass of the outgoing leptons. 



For other phenomenological consequences of non-unitarity in the lepton mixing matrix, see Ref. [21] , 
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The expression for Rp is finally given by 



Ek Fk2 G k2 



Rp 



with 
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(13) 



The result of Eq. (|12p has a straightforward interpretation: i™ represents the impact of new 
interactions (absent in the SM), whereas G %3 encodes the mass-dependent factors. Notice however 
that all states (charged and neutral fermions) do not necessarily contribute to Rp: this can be 
seen from inspection of G 13 , which must be a positive definite quantity. In particular, we denote 
by the N th heaviest neutrino mass eigenstate which is kinematically allowed. 



The SM result can be easily recovered from Eq. (I12j) . in the limit m 
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to which small electromagnetic corrections (accounting for internal bremsstrahlung and structure- 
dependent effects) should be added [T]. 

The general expression for Arp now reads 



Arp 



\{ m P 



m 



2)2 v^-N, 



(0 



max jpml /~mil 



m\{mp 



mf) 2 



Z^n=l 



1. 



pm2 Qn2 



(15) 



Thus, depending on the masses of the new states (and their hierarchy) and most importantly, 
on their mixings to the active neutrinos, Arp can considerably deviate from zero. In order to 
illustrate this, we consider two regimes: in the first (A), all sterile neutrinos are lighter than the 
decaying meson, but heavier than the active neutrino states, i.e. mff tlvc m Vs < mp; in the 
second (B), all v s are heavier than mp. Notice that in case (A), all the mass eigenstates can 
be kinematically available and one should sum over all 3 + N s states; furthermore there is an 
enhancement to Arp arising from phase space factors, see Eq. (I13p . 

We further emphasise that scenarios (A) and (B) are in general experimentally indistinguish- 
able concerning lepton flavour universality, the only exception corresponding to a very particular 
regime where the sterile neutrinos are very close in mass to the decaying pseudoscalar mesor0. 



3 Constraints on sterile neutrinos 

There are strong experimental and observational bounds on the mass regimes and on the size of 
the active-sterile mixings that must be satisfied. Firstly, it is clear that present data on neutrino 
masses and mixings [22] should be accounted for. Secondly there are robust laboratory bounds 
from direct sterile neutrinos searches [15|,l23j . since the latter can be produced in meson decays 
such as 7r — > fJ-^iy, with rates dependent on their mixing with the active neutrinos. Negative 
searches for monochromatic lines in the muon spectrum can be translated into bounds for m Us —9i a 

*In such a situation, the resulting charged lepton would either be less energetic and not pass the experimental 
kinematical cuts [17], or then have a clearly reduced momentum. 
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combinations, where 0i a parametrizes the active-sterile mixing. The non-unitarity of the leptonic 
mixing matrix is also subject to constraints: the rates for leptonic and hadronic processes with 
final state neutrinos depend on ]T\ \ Ul l \ 2 , where (as mentioned above) the sum extends over all 
neutrino states kinematically accessible (i = 1, . . . , -/V max ), and thus constrain the departure from 
the unitarity limit \ Ul l \ 2 = 1. Bounds on the non-unitarity parameter 57 (Eq. (J7|)), were derived 
using Non-Standard Interactions |24| : although not relevant in case (A), these bounds will be taken 
into account when evaluating scenario (B). 

Unless the active-sterile mixings are negligible, the modified Wlv vertex may also contribute to 
lepton flavour violation (LFV) processes^, with potentially large rates, fi — > ej decays, searched 
for by the MEG experiment [26], are the most stringent ones0- the rate induced by sterile neutrinos 
must satisfy [281125] 

BR(^ e 7 ) = "ffirt * 2A >< 10 ^ 12 > ( 16 ) 

where = J2i U^U^G^-^-), with G 7 the loop function and U v the mixing matrix defined 

in Eq. (|6j). Similarly, any change in the Wlv vertex will also affect other leptonic meson decays, in 
particular B — > iv\ the following bounds were enforced in the analysis: BR(i? — > ev) < 9.8 x 10 -7 , 
BR(B -»• pv) < I0~ 6 and BR(B -»• tv) = (1.65 ± 0.34) x 10" 4 [30]. 

Important constraints can also be derived from LHC Higgs searches [31] and electroweak 
precision data [32]. LHC data on Higgs decays already provides some important bounds when the 
sterile states are slightly below 125 GeV (due to the potential -ff-decays to left- and right-handed 
neutrinos). The active-sterile mixings can introduce small deviations to the electroweak fits, which 
allows to constrain them. An effective approach was applied in [32j , assuming very heavy sterile 
neutrinos, and thus these bounds will only be applied in scenario (B). 

Under the assumption of a standard cosmology, the most constraining bounds on sterile neu- 
trinos stem from a wide variety of cosmological observations [15j[33] . Using Large Scale Structure 
(LSS) data, one can also set relevant bounds on very light sterile neutrinos (m Ug < 100 eV), since 
if such light states constitute a non- negligible fraction of the dark matter of the Universe, then 
structure formation is affected. Active-sterile mixing also induces radiative decays Vi — > Vjj, well 
constrained by cosmic X-ray searches. Lyman-a limits, the existence of additional degrees of free- 
dom at the epoch of Big Bang Nucleosynthesis, and Cosmic Microwave Background (CMB) data, 
also allow to set additional bounds in the m Ug — 9i a plane. However, all the above cosmological 
bounds can be evaded if a non-standard cosmology is considered. In fact, the above cosmological 
constraints disappear in scenarios with a low reheating temperature [3l] . In our numerical analysis 
we will allow for the violation of the latter bounds, explicitly stating it. 

4 Arx in the inverse seesaw model 

Although the generic idea explored in this work applies to any model where the active neutrinos 
have sizeable mixings with some additional singlet states, we consider the case of the Inverse 

§ LFV is typically dipole dominated when the sterile neutrinos are light (m Vs < 300 GeV), so that fj, — > cy is the 
most constraining LFV observable. For heavier sterile neutrinos, other (model-dependent) contributions beyond 
the dipole might be more relevant |25| . 

^Recently, it has been also noticed that in the framework of a low-scale type I seesaw, the expected future 
sensitivity of /1 — e conversion experiments can also play a relevant role in detecting or constraining sterile neutrino 
scenarios in the 2 GeV - 1000 TeV mass range |27] . 
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Seesaw [20] to illustrate the potential of a model with sterile neutrinos regarding tree-level con- 
tributions to light meson decays. As mentioned before, there are other possibilities [181119] . 



4.1 The inverse seesaw model 

In the ISS, the SM particle content is extended by tir generations of right-handed (RH) neutrinos 
vr and nx generations of singlet fermions X with lepton number L = —1 and L = +1, respec- 
tively [20] (such that nR + nx = N s ). Even if deviation from unitarity can occur for different 
values of tir and nx, here we will consider the case ur = nx = 3. The lagrangian is given by 

Ass = C SM + Y^D Ri LjH + M Rij v Ri Xj + -n Xij X?Xj + h.c. (17) 

where i, j = 1, 2, 3 are generation indices and H = io^iT*. Notice that the present lepton number 
assignment, together with L = +1 for the SM lepton doublet, implies that the "Dirac"-type right- 
handed neutrino mass term Mr., conserves lepton number, while the "Majorana" mass term fi Xij 
violates it by two units. 

The non-trivial structure of the neutrino Yukawa couplings Y v implies that the left-handed 
neutrinos mix with the RH ones after electroweak symmetry breaking. In the {vl,v r ,X} basis, 
one has the following symmetric (9 x 9) mass matrix M, 

I m T D \ 
M = \ m D Mr . (18) 
V Ml a x ) 

Here ran = -j=Y v v, with v the vacuum expectation value of the SM Higgs boson. Assuming 
fix <C nri£> <C Mr, the diagonalization of A4 leads to an effective Majorana mass matrix for the 
active (light) neutrinos [35] . 

ra v ~ m^M^" 1 fixM^ran , (19) 

whereas the remaining 6 sterile states have masses approximately given by M v ~ Mr. 

In what follows, and without loss of generality, we work in a basis where Mr is a diagonal 
matrix (as are the charged lepton Yukawa couplings). Y v can be written using a modified Casas- 
Ibarra parametrisation [36] (thus automatically complying with light neutrino data), 

/2 

Y u = ^-V^ VM R y^K £4mns . ( 20 ) 

where \/ra v is a diagonal matrix containing the square roots of the three eigenvalues of m u 
(cf. Eq. (|19p ); likewise \/m~ is a (diagonal) matrix with the square roots of the eigenvalues 
of M = Mru^M]^. V diagonalizes M as VMV T = M, and R is a 3 x 3 complex orthogonal 
matrix, parametrized by 3 complex angles, encoding the remaining degrees of freedom. 

The distinctive feature of the ISS is that the additional fix parameter allows to accommodate 
the smallness of the active neutrino masses m u for a low seesaw scale, but with natural Yukawa 
couplings (Y u ~ 0(1)). As a consequence, one can have sizeable mixings between the active 
neutrinos and the additional sterile states. This is in contrast to the canonical type-I seesaw, 
where 0(1) Yukawa couplings require Mr ~ 10 15 GeV, thus leading to truly negligible active- 
sterile mixings. 

The nine neutrino mass eigenstates enter the leptonic charged current through their left-handed 
component (see Eq. ([6]), with i = 1, . . . , 9, j = 1, . . . , 3). The unitary leptonic mixing matrix U v 
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is now denned as C/J \M.U U = diag(mj). Notice however that only the rectangular 3x9 sub-matrix 
(first three columns of U u ) appears in Eq. ([6]) due to the gauge-singlet nature of vr and X. 

In the ISS limit (fix <C m/j <C Mr), and following [37], one can expand the neutrino mass 
matrix in powers of e = thd M^ , block-diagonalizing it at leading order in e, thus easily obtaining 
U v , and the relevant active-sterile neutrino mixing angle. 



4.2 Numerical evaluation of Arx in the inverse seesaw model 



We numerically evaluate the contributions to Rk in the framework of the ISS and address the two 
scenarios discussed before, which can be translated in terms of ranges for the (random) entries of 
the Mr matrix: scenario (A) (m Us < mp) - Mr^ € [0.1,200] MeV; scenario (B) (m Us > mp) - 
Mpi € [1, 10 6 ] GeV. The entries of fix have also been randomly varied in the [0.01 eV, 1 MeV] 
range for both cases. 

The adapted Casas-Ibarra parametrisation for Y u , Eq. (|20p . ensures that neutrino oscillation 
data is satisfied (we use the best-fit values of the global analysis of |22j . and set the CP violating 
phases of C/pmns to zero). The R matrix angles are taken to be real (thus no contributions 
to lepton electric dipole moments are expected), and randomly varied in the range Q{ £ [0, 2ir]. 
Although we do not discuss it here, we have verified that similar Arjc contributions are found 
when considering the more general complex R matrix case. 

In Figs. [U we collect our results for Ar^ in scenarios (A) - left panel - and (B) - right panel, 
as a function of fj, which parametrizes the departure from unitarity of the active neutrino mixing 
sub-matrix C/pmnSj fj = 1 — |Det(C/pMNs)|- Although the cosmological constraints are not always 
satisfied, we stress that all points displayed comply with the different experimental and laboratory 
bounds discussed before. For the case of scenario (A), one can have very large contributions to 
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Figure 1: Contributions to Ar^- in the inverse seesaw as a function of fj = 1 — |Det(C/pMNs)l : 
scenarios A (left) and B (right). The upper (lower) dashed line denotes the current experimental 
limit (expected sensitivity). On the right panel, red points denote cases where Y v > IO" 2 . All 
points comply with experimental and laboratory constraints. Points in (B) are also in agreement 
with cosmological bounds, while those in (A) require considering a non-standard cosmology. 

Rk, which can even reach values Ar^ ~ 0(1) (hi some specific cases we find Arx as large as 
~ 100). The hierarchy of the sterile neutrino spectrum in case (A) is such that one can indeed have 
a significant amount of LFU violation, while still avoiding non-unitarity bounds. Although this 
scenario would in principle allow to produce sterile neutrinos in light meson decays, the smallness 
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of the associated Y v (< O(10~ 4 )), together with the loop function suppression (G 7 ), precludes the 
observation of LFV processes, even those with very good associated experimental sensitivity, as 
is the case of fi — > cy. The strong constraints from CMB and X-rays would exclude scenario (A); 
in order to render it viable, one would require a non-standard cosmology. 

Despite the fact that in case (B) the hierarchy of the sterile states is such that non-unitarity 
bounds become very stringent (since the sterile neutrinos are not kinematically viable meson decay 
final states), sizeable LFU violation is also possible, with deviations from the SM predictions again 
as large as Arx ~ 0(1)- Contrary to case (A), whose results could also arise in other frameworks 
with light sterile neutrinos, the large deviations in (B) typically occur when all the singlet states 
are considerably heavier than the decaying meson, and reflect specific features of the ISS. As can 
be inferred from Eq. (|19p . in the inverse seesaw framework, one has m u ~ (Y u v/Mr) 2 fix] hence, 
for "low" (when compared to, for instance, the type I seesaw scale) Mr, light neutrino data can 
still be accommodated with large Yukawa couplings, Y u ~ few x 10 _1 . As a consequence, large 
active-sterile mixings can occur, thus leading to an enhancement of Rk- Even if in this case one 
cannot produce sterile states in meson decays, the large Y v open the possibility of having larger 
contributions to LFV observables so that, for example, BR(/x — > ej) can be within MEG reach in 
this case. 

Although we do not explicitly display it here, the prospects for Arv are similar: in the same 
framework, one could have Ar n ~ O(Arx), and thus Ar n ~ 0(1) in both scenarios. Depending 
on the singlet spectrum, these observables can also be strongly correlated: if all the sterile states 
are either lighter than the pion (as it is the case of scenario (A)) or then heavier than the kaon, 
one finds Ar n ~ Ar^. The latter possibilities are a feature of the ISS mechanism (not possible in 
the unconstrained MSSM, for example) and are expected to be present in other low-scale seesaw 
models that allow for large active-sterile mixing angles. 

5 Concluding remarks 

The existence of sterile neutrinos can potentially lead to a significant violation of lepton flavour 
universality at tree-level in light meson decays. As shown in this study, provided that the active- 
sterile mixings are sufficiently large, the modified Wiv interaction can lead to large contributions 
to lepton flavour universality observables, with measurable deviations from the standard model 
expectations, well within experimental sensitivity. This mechanism might take place in a number 
of frameworks, the exact contributions for a given observable being model-dependent. 

As an illustrative (numerical) example, we have evaluated the contributions to Rk in the 
inverse seesaw extension of the SM - a truly minimal extension of the SM - , for distinct hierarchies 
of the sterile states. In particular, we have studied the impact of non-unitarity in a low mass regime 
for the additional singlets, an inverse seesaw mass regime considerably lower than what had been 
previously addressed [211124] . Recent studies [38] have proposed a search of the monochromatic 
peak in the next generation of high intensity experiments, yielding both the mass and mixing 
angles for sterile neutrinos with masses in the range 3 MeV < m Va < 414 MeV. 

Our analysis reveals that very large deviations from the SM predictions can be found (Ar^ ~ 
0(1)) - or even larger, well within reach of the NA62 experiment at CERN. This is in clear 
contrast with other models of new physics (for example unconstrained SUSY models, where one 
typically has Atk ^ 0(1O~ 3 )). We further notice that these large deviations are a generic and 
non fine-tuned feature of this model. It is worth emphasising that, in view of the potentially large 
new contributions to these observables, such an analysis of LFU violation in light meson decays 
actually allows to set bounds on the amount of unitarity violation (parametrized by 77). 
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Interestingly, in this framework, both Arj< and Ar n are strongly correlated in the case where 
all the sterile states are lighter than the pion or heavier than the kaon. 

The impact of this mechanism is not restricted to light meson decays: there are currently 
some hints of lepton flavour universality violation in heavy mesons, with deviations already found 
in observables such as T(B~ — > tu)/T(B° — > tt + £~v). We expect significant contributions to 
.B-meson observables [39], which have very promising experimental perspectives. 
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